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ABSTRACT 


A two-phase air-water propulsion system was studied to establish the 
scaling laws and to determine the efficiency. A calculational method 
which separates frictional, residuary and air flow associated forces was 
postulated and physical models were built and tested in a towing tank to 
provide data for verification of the postulate. Because of the poor qual- 
ity of the data, the strongest statement that can be made at this time is 
that the postulated method of separation of forces holds promise and that 
the system can be expected to have an efficiency (exclusive of power plant 
and piping, etc.) of about 20%. 


Some of the special features and potential .areas of application of 
the two-phase propulsor are discussed. The experimental methodology was 
examined in light of the results achieved by its use and it was concluded 
that the experimental procedure needed improvement. Suggestions for fur- 
ther work in the field are given. 
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In recent years several novel ship propulsion systems have been 
proposed, tested and some have been utilized in full scale. One of the 
propulsion systems that has been proposed recently, henceforth referred 
to as a two-phase propulsion system, is the subject of study of this 
thesis. 

Mr. William Fulton has applied for patent rights for this unique 
propulsion system which consists of air under pressure being forced 
through an orifice, located approximately amidships at the bottom of 
ths ship, into a channel which conducts the air to the stern of the 
ship. The bottom and the sides of the aft part of the ship and the 
water beneath the ship form the sides of the channel. The channel is on 
an angle so that as the air-water mixture moves aft it also rises toward 
the surface of the water. 

Professor 5. Schuster was the first to study the theory of the two- 
phase propulsion system and published his results in a Professor 
Kostilainen has done considerable experimental work on this type of pro- 
pulsion system and has repcrted his results in Several papers. 
The results of Kostilainen's work are restricted in application but the 
work has accomplished the 515: 

(1) The dimensional analysis of the problem and the establishment 
Of some useful dimensionless ratios. 

(2) The calculation of the momentum and energy balance equations 


which identify some of the functional relationships which govern the 


phenomenon under investigation. 





= He 

(3) The measurement of the air-flov rate and the apparent resist- 
ance of the model with the propulsor operating for two geosims which had a 
scale factor of four. (The waterline length of the small model was 1.6 
meters.) The apparent resistance differs from the total resistance by 
the amount of thrust developed by the propulsor. 

(4) The generation of data and presentation of some of the theory 
of bubble formation and spreading in the channel. 

The major restriction of Kostilainen's work is that although models 
of two different sizes had been studied, no scaling law was postulated 
by Kostilainen. In addition the results published do not allow other re- 
searchers to postulate such a scaling law and be in a position to verify 
it. In view of the fact that what is of importance is the behavior of 
the full sized ship, it was concluded that a study that would permit the 
establishment of such a scaling law was fuliy justified. 

The proposed study involves: 

(1) The formation of. a mathematical model that interrelates the 
critical parameters that describe the phenomenon under investigation and 
permits the prediction of the full scale ship. 

(2) Model testing is to be used to verify the proposed mathematical 
model. 

(3) The comparison of the performance of the proposed propulsion 
system with the known performance of other propulsion systems. 

should be noted that attempts to calculate thrust from first 
principles of hydrodynamics were not successful. Therefore it was de- 
cided that an experimental approach would be used. 


Mr. Fulton had a rough, three-foot model which he had used as the 
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basis for his patent application and which he was willing to provide for 
use in this thesis work. His model was faired and slightly modified and 
a geosim (Model A) with a 1.5 scale factor was built and the two models 
were tested without air in the M.I.T. towing tank. (For details of 
Model A see Ship Model, Chapter IV, Section A.) The tests proved that 
the original model which was operating at a Reynold's number of between 
1075 and 8 x 1057 could not be driven into a turbulent flov and the result- 
ing drag measurements could not be compared to those of the larger model. 
Therefore an even larger model (Model B) was built, again on a scale fac- 
tor of 1.5. The two models, A and B, used in the final comparison tests, 
were therefore 4.5 feet and 6.75 feet in length respectively. 

Figure 1 shows the geometric configuration studied in this work. 

The model is bottom up in the pictures. Air under pressure enters the 

air chamber and passes out the astern opening into the channel. As it 
passes up the channel toward the stern it imparts momentum to the adjacent 
water and thus creates a reaction force on the ramp, thus propelling the 
model. 

Observations made through the glass wall of the towing tank indicate 
hət two modes of operation occur. One, at low air flow, is characterized 
by bubbles passing up the channel. The second, at higher air-flow rates, 
is characterized by the establishment of waves traveling up the channel. 
Ub course, there is a range of air flows which displays characteristics of 
both these modes. If the air-flow rate is increased to a high enough 
value some of the air will start 25777 from the air chamber and rise 
up Outside the ship and provide no propulsion force. 


Chapter II gives a more detailed study of the theoretical and 











Eq 
experimental results which have been published to date. Chapter III 
presents the theory used in organizing this work. Chapter IV describes 

the equipment used including the model, the tank, the force detectors 

and the pressure detectors. In addition it details the way in which 

each piece of equipment was used and the analytic procedures used to ob- 
tain the results. Chapter V presents the detailed results of the work 

in graphical or tabular form. Chapter VI draws upon the results of the 
preceding section and presents the conclusions that were reached concern- 
ing the results themselves and the experimental procedure used in obtain- 
ing them. Chapter VII gives some potential applications of the two-phase 
propulsion system and gives a listing of some of the engineering considera- 
tions and problems which might arise from use of this propulsor. This 
7:5: also provides someone interested in further work on the two-phase 
propulsor with some ideas about what has not yet been investigated but 
deserves further consideration. 

Appendix A is provided for reader convenience and simply details 
Kostilainen's momentum balance and energy balance calculations. Appendix 
B shows the details of calculating Cy > Cro and n (the principle parameters 
used in describing the phenomenon under investigation) for one cperating 
condition. 

Tables in this thesis provide the averaged data from the experiments 
so that it will be available to future researchers that wish to adopt a 
different analysis approach to the problem than the one followed by the 


Embhor. 
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dn STATE OT THE AT 


A. Theoretical 

As noted in the Introduction, Professor Kostilainen of the Finland 
İnstitute of Technology has been studying the tvo-phase propulsion system. 
He has completed a dimensional analysis and solved a momentum balance and 
an energy balance on an idealized model of the system. The details of his 
‘calculations can be found in Appendix A, but the results are summarized 
below. The notation in the equations below is not Kostilainen's original 
notation but rather that of the author. The change in notation was 
prompted by the fact that Kostilainen has used some symbols which have 
very specific meaning to naval architects and stand for different quanti- 
ties. For example he calls B the breadth of the channel, but B is most 
commonly used to denote ship's bean. 

Kostilainen's dimensional analysis results in the following three 


equations. 


“TS “ pgBch2 


EN s: C. = Thrust coefficient 


— 
" 


Thrust created by the tvo-phase propulsion system, parallel 
to the ship"s fore and aft axis 


Mass density of vater 


© 
11 


Acceleration due to gravity 


09 
H 





1 
P 1لا‎ 65 of the channel 


h = Depth of the orifices 


$, - Air flow rate coefficient 
Ve = Volumetric air flow rate 
D = Froude number (modified) 


UNE Veleros theMiiquia at the orifices 
Unfortunately the results of his momentum balance include quantities 
which cannot be measured, such as the amount of water that gets acceler- 
ated and the final velocity of the water leaving the ramp, and are there- 
fore not very useful to an experimenter. He does however show that the 
propulsor efficiency, n, is 


TUL 
pVcgh 





n = 





007 
D. Experimental 

| Published information from Kostilainen to date describes two geosim 
models with a scale factor of 4 and the results of his self propelled 
tests. Because he decided not to measure the thrust of the two-phase 
propulsor he wee obliged to find another way of presenting his data, and 
177 -reTfore defined an "apparent thrust," TA, which is the force re- 
quired to tow the model at a prespecified speed without any air flow, 
henceforth referred to as the bare hull drag, minus the force required 
"tov" the model at the same speed with the propulsor in operation. It 
should be pointed out that the apparent thrust includes not only the real 
propulsor thrust but also any changes in the actual hull resistance which 
result from changes in the flow pattern and pressure distributions caused 
by the presence of the air in the channel. The apparent thrust coef- 


ent then is 


Ta 
PTA = 5 
PgBch 


Kostilainen found, using the above dimensionless coefficient, that 
for low air flow rates, i.e., 9 > .2, The bollard pull tests, i.e.. Tess. 
with model velocity of zero, the results were independent of model size 
and at higher flow rates the bollard pull test results showed weak depend- 
ence on model size. The dependence on model size increased with Froude 
number and at the self propelled point was large.enough so that correla- 
tion was almost nonexistent. 

Some very important results of Kostilainen's work are that the model 
ШОС Ту, Ј, 15 а fair approximation to Ü,, the liquid velocity at the 


orifices, and that the thrust is mainly dependent upon the Froude number 











-22- 


ఏ హక 


Vhen ship propulsion is being considered at least three quantities 
are of interest: thrust from the propulsor, its efficiency and resistance 
Of the hull. Hull resistance for a normal displacement vessel is knovn 
to be fairly vell predictable by Froude"s hypothesis and it is mostly a 
function of the Froude number and the Reynolds number. 

For a tvo-phase propulsor driven ship it is postulated that the hull 
resistance is a function of the dimensionless air flov rate, the Froude 


number, FP, and the Reynolds number, R, i.e., 


Ue xcu 





Cm = bpSU?f —— — 
T L?VEL VEL V 
vhere 
Co = total resistance coefficient 


p = density 
S = wetted surface of the ship 
U = speed of the ship 
77 ee metric air flow rate 
g = gravitational constant 


య ఆ or ship 


v = viscosity of water 


It is additionally hypothesized that: 





Cr V 

U ЫШ G U 
- 0 с + C P3— 5?) -— 2 2 
15050 dE e |== vgl 


where 


Cp = residuary resistance coefficient 





o 





805 ——-. 


200 


Cp frictional resistance coefficient 


a rli resistance ceefficıent 


Co 
$ is now defined as 


Ус 
LÖVEL 





Cp Plus Cr is interpreted as the bare hull resistance, i.e., the 
resistance with no air flowing in the channel. Having made this final 
postulate, the resistance equation can be interpreted as saying that the 
bare hull resistance, Cry > plus the air flow resistance equals the total 
resistance, i.e., 

27: “ән 

If all of the above is correct, then a knowledge of Cg as a function 
Of the dimensionless air flov rate and the Froude number vill allov pre- 
diction of the full scale ship resistance from the model resistance. 


Knowledge of C,, alone is not sufficient information to do a powering 


T 
calculation for the tvo-phase propulsion system. Two options will be pre- 
sented in this thesis for obtaining the power requirements of a full sized 
“hip. 

The first option is to try Kostilainen's nondimensionalizing coef- 
ficients on the thrust (not the apparent thrust as he did). Therefore a 
new thrust coefficient, CTS» is defined as 

CTS = ү 
Since it does not contain the change in hull resistance caused by air in 


the channel as Cra did in Kostilainen's work, 1t may correlate better than 


Cr, did. 





m. 

Under this first option then, Co and Crg are presumed to be known 
for the model and they are of course both functions of q and F. Along 
with knowledge of Cy and Стс it is assumed that the Cpy has also been 
Measured for the model so that Cr and Cp are known. To obtain the power- 
ing requirements for a full sized ship with a known scaling factor and 
at a specified speed the following procedure could be used. 

The scaling factor and the speed can be used to calculate R and F 
Men can be used to obtain Cp and Cr for full scale, thus providing Cpy 
full scale. Entering the С̧у curves vith the given F provides Cə .values 
for various values of q. These values of Cg are added to Cpy to give Ст 
from which the total hull resistance can be obtained and then plotted as 
a function of $. Now Cyc can be obtained from its curves as a function 
Of $, and, since pgB..h” for the full sized ship is known, the thrust can 
be calculated and plotted on the same graph as the resistance was plotted. 
Where the curves cross the thrust equals the resistance and the ship is 
self propelled. The crossover point also determines $ and therefore the 
input power requirement can be calculated since 

Power Input = Vopeh. 

The second option is that of recognizing that the efficiency, n, is 
a ME ~: ~: number and may be a useful tool in extrapolating from 
the model to full scale. Since the power ote and output can be calcu- 
lated from measurable quantities on the model, it may be possible to es- 


tablish a family of curves of n as a function of $ and F. 


5 uv 
Vepgh 





So for this option it is presumed that Co and n are known as a function of 





05 
ç and F. To obtain the powering requirements for a full sized ship with 
a known scaling factor and at a specific speed the following procedure 
could be used. 

The total hull resistance can be obtained as it was in the preceding 
option and plotted against 6. The thrust can be calculated for the full 
scale by entering the family of n curves with the specified F and thereby 
obtaining n as a function of q. Utilizing that knowledge and the equation 
above for n it is possible to calculate T for specific values of ¢ and 
then plot them on the resistance plot. At the crossover point the resist- 
ance equals the thrust and the ship would be self propelled. © would then 
be known at the self propelled point and the input power requirement cal- 


culated. 
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IV. EXPERIMENTAL PROCEDURE 


A. Equipment 

Two geosim models were built on a scale factor of 1.5 having L = 4.5 
feet and 6.75 feet. Figure 2 is a lines drawing of the geometric 
configuration used. Table 1 lists the major dimensional character- 
istics of the two models. A is defined as the vertical height of the air 
chamber. 46 is the inner diameter of the air supply pipe to the air cham- 
ber. The models were made of six pound per cubic foot urethane resin 
insulation which was cut and molded to shape and then sealed with three 
coats of fiberglas finishinz resin. Each coat of resin was sanded to a 
smooth finish. Two coats of spray epoxy paint were then used to finish 
the surface. Turbulence pins were glued to the models at both stations 
2 and 4 on a 1/2" and 3/4" spacing between pins for the small and large 
model respectively. The force detector mounting platforms were located 
with their midpoints at station 8 on the waterline. Lead ballast was 
used to give each model its proper displacement and zero trim. 

The models were towed in the Massachusetts Institute of Technology 
an Model Towing Tank. The tank width is 102 inches and the water 
depth vas 42 inches. The carriage travel had a programed acceleraticn 
followed by 150 foot travel at uniform speed and then a programed deceler- 
ation. Speeds are accurate to within 1%. Over the period of time the 
tests were conducted the fresh water in the tank did not vary over one 
degree from 78° F. All calculations have been based on that temperature. 

Two different force detectors, a torque meter and a force block, 
were used because the torque meter could not tolerate the large accelera- 


tion forces that were experienced during the acceleration and deceleration 
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MODEL PARAMETERS 


Description 


Length 

Breadth 

Breadth of channel 

Wetted surface area 

Vertical height of air chamber 
Angle of ramp to water surface 


Inner diameter of the air 
Supply pipe 


Ef: 16 Coefficient 
Depth of orifice 


Displacement 


Model A 


hoo WT. 


POS 


1007511 


SS 기 


Tone 


12.19 


చ 
. 996 
5 625 I 


75.3# 


Model B 


6.7/5 ee 
15255 im? 
ages 
17.50 sq.ft. 
Tn 


12:14: 


To ири 
. 996 
SHH ane 


254.3# 





-29- 
of the carriage when the large model was being towed, but it gave greater 
accuracy for measuring the smaller forces associated with the smaller 
model. Both meters were calibrated several times during their use and 
each was determined to be linear to within 1% to Forces onunu ui 

3 1bf for the torque meter and force block respectively. 

The torque meter was a variable resistance balanced bridge type 
while the force block was a variable reluctance detector. Both were con- 
nected to a transducer with a D.C. voltage output which was battery biased 
with a 5 vde battery. The biasing vas used to keep oscillations.in the 
D.C. level from passing through zero thereby causing malfunction of the 
voltage to frequency converter vhich vas used to feed a pulse counter 
that vas automatically started by the carriage after it had reached uni- 
form speed. 

The model was further instrumented with four pressure taps feeding 
pressure transducers. The taps, P1, P5, P4, were located at three equal- 
ly spaced points on the channel ramp so that P, was nearest the air cham- 
ber and P4 was nearest the stern and they divided the ramp into four equal 
sections. The fourth tap, Py, was installed in the air supply piping 
st above the air chamber. The pressure transducers had a 1% accuracy 
over the range 0 to .5 psi. They were of the variable reluctance balanced 
bridge type used with carrier preamplifiers which provided resistance and 
capacitance balance plus rectification of the output signal. The output 
was recorded on a multichannel chart recorder. The pressure recording 
system was calibrated several times during the course of the experimenta- 
tion using a manometer and a static head of up to ten inches of water. 


Air was supplied to the model from 130 psi air system on the carriage. 
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A pressure reducer vas used to control the pressure at the inlet to the 
rotometer used for measuring the air flov. The air flov rate vas set 
by means of a needle valve. The rotometers vere of the gravitational 
floating ball type and had maximum flow rates of 92 cfh and 1482 cfh. 
The air supply could provide only about 700 cfh without significant loss 


of pressure and therefore loss of accuracy of measurement. 





oe 
Bi Operational Procedure 

The operational procedure, i.e., those functions which were per- 
formed to obtain the data, were for all practical purposes the same for 
both models. 

The first requirements were those of calibration. The force detector, 
after being connected to the transducer, the battery, the voltage to fre- 
quency converter and the counter, was calibrated by using a one-pound 
weight to supply the force. Both force detectors were also checked for 
linearity by using various weights over the ranges of interest. „Both de- 
tectors were linear to within +1%. 

The pressure transducers also required calibration. They were all 
connected to a water manometer anc their output measured on the recorder, 
after resistive and capacitive balancing had been accomplished on the 
electrical system. The manometer readings were plotted against the re- 
corder readings for various pressures from one to ten inches of water to 
check linearity and for accurate calibration. The recorded outputs were 
linear to within +5%. 

Each model was attached to its appropriate force detector, which was 
attached to the carriage. The model had freedom of motion in both pitch 
and — but was restrained in yaw by the force detector. After the 
model was attached to the carriage the air supply flexible piping and the 
pressure taps flexible connections were made and sealed using rubber seal- 
ant. The transducers were again balanced both resistively and capacitive- 
ly and the static pressure was nulled out by using bias on the preamp so 
that the detector measured the deviation from the static condition. 


Bare hull resistance runs, i.e., without any air flowing, were then 
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made and the force, speed and pressure distribution measured. Each run 
consisted of running the model the length of the tank at uniform velocity 
(exclusive of acceleration end decelerations at the beginning and end of 
each run). Several runs were made for each speed. 

Bollard pull tests were conducted with the carriage stopped, i.e., 
zero velocity, and the air flowing at some desired rate. The air flow 
rate was measured with the rotometer and the force and pressure distribu- 
tiontion were measured as in the bare hull resistance runs. 

If the air flow rate changed over the run by more than 2% the run 
was performed over again. The bollard pull test was conducted several 
times on each model so that the average results could te used with confi- 
dence of repeatability. 

Finally, runs which henceforth will be called the self-propelled 
model tests, were conducted near the self-propelled point and the force, 
velocity and air flov rate were measured and the pressure distribution 
was recorded. The procedure was to start the air flowing and adjust it 
to the desired rate with the model stopped. The model was then acceler- 
ated to speed and the force on the detector measured and the pressure 
ШЕ ution recorded. After the model was stopped the air flow rate was 
again measured and if the rate had changed by more than 2% the run was 
conducted again. Several runs were made at each velocity-air flow rate 
combination. 

After the self-propelled tests had been completed at a given velocity, 


the bare hull resistance was measured again at that velocity. 





5: 
. Analytical Procedure 
(1) Bare Hull Resistance Coefficients. 
The first calculations to be made were those which established 
that the models operating without air obeyed Froude's hypothesis in spite 


of the presence of the ramp in the aft section of the ship. 


DT 
LSU? 


- RE Cp = 

where Dp = total resistance (measured). 
Reynolds number anc Froude number меге calculated for both models at 78% F. 
Knowledge of Reynolds number provided the required input for calculating 


( 
Cr accordance with the formula for the ITTC line 


Ор = 0.075/(l0g,¿R - 2) 


The average of several runs at the same speed was used to increase 
the accuracy of Dp and Cp was then easily calculated. Cp was simply then 
the difference between Co and Cr. “For the model van nu ıı moe n 으 
be used after the above calculation vas complete but for the longer model, 
L = 6.75 ft., a correction had to be made because the model vas large 
enough so that the model-wall and bottom interaction force was large 
enough to affect the resistance coefficient. This interaction force, 
after being nondimensionalized by division by LoSU?, is called the tank 
blockage coefficient, Cj. | 

5 has presented a vay for calculating Cg based on the tank 
dimensions, the model displacement and dimensions and the Reynolds number. 
These numbers have been presented in Table 3. Their inclusion in the 
above calculation then resulted in values of Cp for the large model (L = 


6.75 ft.). Cp is plotted and discussed in the next chapter. 
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(2) Bollard Pull Tests 

The next calculation that needed to be made was the bollard 
pull calculation which is actually two separate calculations. One part 
of the calculation is very trivial and simply required that the forces 
as measured by the force detector be averaged over the runs that were 
made at a given air flow rate. This force could then be nondimensional- 
ized by division by the factor రం if desired to yield Cro. 

The second part of the bollard pull test was more complicated to 
perform but the idea behind it was quite simple. The pressure distribu- 
tion on. the ramp and the front of the air chamber, called the step, was 
to be integrated to give the total force created by the two-phase air- 
water mixture. | 

The pressures at the center of the channel had been recorded on 
charts and therefore had to be first converted from chart dimensions to 
pressures using the various calibration factors that had been measured. 
The pressures were then averaged over all the runs which had the same 
air flow rate. It was assumed that the pressure distribution was not a 
function of the athwartships position. Plots of these pressures are in 
the results section, Figures 6 and 7. The area under each curve is 
proportional to the force on the ramp and therefore each curve was inte- 
grated using Simpson's first NE The force on the ramp is equal to 
the sine cf the ramp angle times the width of the channel times the 
integral of the pressure distribution curve (which has units of pressure 
times distance). 


The above integrated force, F,, does not account for the force cre- 


R 
ated by the over pressure in the air chamber acting against the step. 


This force is called the step correction force, Step Corr. The over 





Sage 
pressure in the air chamber is P, but it does not act against the full 
area of the step, only an area which is equal to the width of the step 
(channel width) times the over pressure in inches. The idea here is 
that the over pressure, Pus in the chamber is on the average equal to 
the vertical height of the layer of air in the chamber since En equal to 
zero implies a static pressure of h - A inches of water. One indication 
that this is a more reasonable approach to the calculation of the step 
correction force than that of using the full step area is that calcula- 
tions made using the full area resulted in calculated forces greater than 
the measured forces. 

By adding the ramp force and the step correction force together 
the total force created by the two-phase 270075” could be obtained, T. 
This force was used in both a dimensional and nondimensionel form. See 


Figures 8 and 9. 


(3) Ramp Drag 
The ramp drag or bare hull ramp drag, Rp; is a name that has 
been given to the final results of the calculation which is described 
ES. Essentially it is the sum of the Step Corr and the integrated 
force on the ramp when the model is moving but no air is flowing. It 
has no real physical meaning and is simply used as a tool to achieve a 


result. 


(4) Thrust Provided by the Two-Phase Propulsion System 
Took ahead it the thrust calculation in Section 175 
disclose that under many velocity-air flov rate combinations the pressures 


are mostly negative and that the result of the integration is negative yet 
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there has been a decrease in the measured force required to move the hull 
at that speed compared to the bare hull runs. The point is that the in- 
tegrated pressure on the ramp is less negative with air than without and 
the net difference is the thrust on the ramp provided by the two-phase 
propulsion system. This thrust could have been calculated by taking the 
pressure distribution with air and subtracting the pressure distribution 
without air from it and then performing an integration procedure like the 
one detailed in the preceding section for the bollard pull tests. How- 
ever, since the integration is linear in pressure and distance the two 
pressure distributions can be done separately and the results subtracted 


to yield the same answer. 


(5) Self-Propelled Model Tests 
Self-propelled model tests are to provide three additional pieces 
Of information: the self-propelled speed as a function of the air flow 
rate, the resistance of the hull vith air floving, and the thrust from the 
two-phase propulsion system. 

The analytic procedure required to get the self-propelled speed as a 
Eon of the air flow rate could have been as simple as plotting the 
counts measured by the counter measuring force against the air flow rate 
and then establishing the constant velocity lines. When the counts passed 
through zero (i.e., at zero force) the model had to have been self propel- 
led at that air flow rate and velocity. To make 01 graphs more meaningful 


¡en 


(D 


counts were changed to forces. For accuracy the results of many runs 
for each velocity-air flow rate combination were averaged and then plotted. 


To present the data in a more easily assimilated form, the self- 





2 
propelled velocities and air flov rate combinations for each model vas 
nondimensionalized to F and $ and plotted together. (See Figure 15.) 
These numbers vill be needed later to locate the self propelled points 
277 Der cürves. 

The Porce, ın measured by the force detector when the model is mov- 
ing with the air flowing when subtracted from the thrust developed by the 
propulsor yields the second desired result, i.e., the resistance of the 
bul, Ву. D and ру can be nondimensionalized as needed. (Note that 
vhen the air is providing more thrust than is needed to overcome the re- 
Excuse of the hull that the measured force T, is positive.) 

The thrust provided by the two-phase propulsion system is to be ob- 
tained by integrating the pressure distribution. The pressures were re- 
corded on chart paper. Calibration numbers were used to change the chart 
dimensions to inches of water and the pressures were then averaged over 
runs mith almost the same air flow rate and velocity. These pressure 
distributions were then handled in the same manner as the bollard pull 
tests were with one exception. To obtain the thrust it is necessary to 
Subtract the ramp drag from the integrated force on the ramp and the step 
correction force. The resultant is the total thrust created by the two- 
phase propulsion system along the ship's fore and aft axis, T. 


The final few calculations to obtain C,(F, $), n(F, 9) and Cos C» $) 


ф 
are simply obtainable from the thrust by arithmetic. 


E "UE a 


BH 
1soSU2 


All of the numbers are known and Cpy has already been calculated as part 


of the procedure of IV.C.1, the calculation of Cp to prove Froude scaling. 
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calculated using 


Bi: TU 


Since the value of all of the terms has been established. Сес 157657 


tained from the thrust using the nondimensionalizing coefficient 1/ogB.h”. 
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Vo EXPERIMENTAL KESULİS 


A. Establishment of Froude Scaling of the Models 

Both models were towed in the tank over as wide a range of speeds as 
the tank equipment and models could permit. Though the steady state 
forces were generally less than two pounds force, the inertial accelera- 
tion forces on the carriage and the force block or the torque meter were 
large enough to cause concern for those pieces of equipment and thereby 
limited the range of speeds available for testing. 

Below a Froude number of 0.08 the scatter in data was unacceptably 
high and therefore the data are not presented. The scatter was caused by 
random transition from laminar to turbulent flow under conditions which 
were such that the frictional contribution to the model's resistance was 
greater than 50%. 

Figure 3 shows the residuary resistance coefficient plotted as a 
function of the Froude number. The total Cp» CF and the blockage correc- 
EX" coefficient if applicable are tabulated in Tables 2 and 3. 

Cr is never the dominant contributor to Cp and therefore the magnitude of 
Every sensitive to the fractional contribution the experimenter pos- 
tulates. The ITTC line was used in that calculation because the history 

of the towing tank indicates that has been the most successful estimation 
OT Cn and given the best correlation with larger models, full sized ships 
and models run at other tanks. 

Without the blockage correction for the 6.75 ft. model the correla- 


tion would be very poor. With it, however, the models are seen to be 


basically similar above a Froude number of about 0.15. Below that value 
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ee 
the 4.5 ft. model is in the transition zone of the fractional contribu- 
tion while the 6.75 ft. model is still well above that zone. This causes 
the calculation procedure to artificially depress the Cp curve of the 
smaller model. 

Also plotted in Figure 3 is the 4.5 ft. model as measured with the 
force block, and this plot is presented as verification that the two 


force measuring devices yield similar results. 
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B. Eollard Pull Tests 

The bollard pull test is of interest for two reasons. It establishes 
the zero velocity value of the thrust. In addition to that needed func- 
tion it can also be used as a source of information concerning errors for 
the forces calculated from pressure-area integrations. The bollard pull 
force and the bollard pull pressure distribution were measured simultane- 
ously and therefore the error introduced by the integration process can be 
evaluated by comparison of the measured and calculated values. 

Figures 4 and ° show the pressure distribution of these two tests. 
Nondimensional plotting of these quantities using the factors that are 
used elsewhere in this thesis does not cause the dissimilarities to dis- 
appear. The most noteworthy dissimilarity is that Py is always the highest 
pressure on the 6.75 ft. model, and this is not the case on the 4.5 ft. 
model. The second most important change is that P4 and P, have changed 
౨౧౧5. Since Pa 15 the closest measurement point to PI TS n s 


EN St FP. is highest with the lowest being Pj. These differences imply 


1 
that the flow patterns of the two models are different. In order to ob- 
tain a little more insight before we try to determine the consequences of 
the above implication the pressures were plotted against their position on 
the ramp. Again the dissimilarities are easily spotted, even though the 
Overall shapes are comparable. 

One of the basic assumptions that has been made in attempting to under- 
stand the two-phase propulsion system is that some form of nondimensional- 
ization can be used to make the hulls "similar" on the basis of force scal- 


ing. It is expected that, even though the flow patterns are dissimilar, 


the forces of interest can be scaled by the use of the proper procedure. 
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Figure 4 


Bollard Pull Test Pressure Distribution as a 


Function of the Air Flow Rate for Model A 
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Bollard Pull Test Pressure Dristribution as a 


Function of the Air Flow Rate for Model B 
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Figuro 6 


Bollard Pull Test Pressurc Distribution 


as a Function of Position for Model A 
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Bollard Pull Test Pressure Distribution 
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ug 
This hopeful expectation is based on the successful application of Froude's 
hypothesis. 

One plausible explanation of the dissimilarities is that frictional 
forces are more predominant in the smaller (L = 4.5 ft.) model. 

Kostilainen has postulated that this is the case and — work seems to 
verify that postulate. The reason for believing this is that the increase 
in dynamic pressure is shifted in the direction of the source on the small 
model, thus indicating that the air-water mixture has been obliged to do 
27 r,, i.e., has lost more energy over the first portion of its travel. 
This also implies that it has proportionally vorked against a larger force. 
This correlates with the fact that on the smaller model the local Reynolds 
number for the position of each of the three probes on the ramp is smaller 
than for the large model which implies that frictional forces will be more 
dominant. 

One conclusion that can be dravn nov is that any nondimensionalizing 
scheme vhich does not treat the frictional forces separately vill probably 
not be successful at scaling the total thrust produced by the two-phase 
propulsion system. 

Figures 8 and 9 show the comparison of the calculated forces to the 
measured forces in each case. The steps which lead from the pressure dis- 
tribution to the integrated force value are tabulated in Tables 4 and 5 
Both integrals yield low values. They are not proportion errors. One ap- 
parently is in error by a fixed additive value, while the other if off by 
a fixed multiplicative value. At least a part of the error is probably 


the result of assuming that the pressure distribution is flat across the 





250- 
Figure 8 


Bollard Pull Test Calculated and Measured Forces 


as a Function of Air Flow Rate for Model A 
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Figure 9 


Bollaxd Pull Test Calculated and Measured Forces 


as a Function of Air Flow Rate for Model B 
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5 
channel. TFrom visual observation it is knovn that some of the flov pat- 
terns in the channel vere such that at the center of the flov pattern the 
water was closer to the ramp than on the edges, i.e., the sides and cor- 
ners were full of air and the center had relatively less. If the system 
is assumed to be in quasi-equilibrium, then the pressure in a bubble will 
be approximately constant and a flow shaped as described above would have 
higher pressures off the center of the ramp than on it. There were flow 
patterns observed which were the opposite of this, i.e., more air in the 
center and this would give a lower pressure off the center line. . The cor- 
relation between these small flow differences and the magnitude of the 
difference between the measured and the calculated forces has not been 
fully determined. 

The two curves in each plot of Figures 8 and 9 cannot be used as a 
source of a "calibration" factors for correcting all the integrated forces 
to their proper value because the pressure distribution is affected by the 
velocity of the model and therefore the correction factor which brings the 
zero velocity integrated force to the correct level is probably not the 
right factor to use at any other velocity. Therefore no attempt will be 
made to modify the integrated forces at other velocities because the modi- 
fication will only bring additional confusion to the results. This lack 
of equality will, however, be alluded to occasionally to explain why some 
results appear as they do. 

By inference from the fact that the inequality exists at zero velocity, 
it is claimed that the calculated forces for the lower velocities are prob- 
ably smaller than the real forces they are to represent. In addition it is 


also claimed that for the smaller model (L = 4.5 ft.) the magnitude of the 
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Es Determination of the Ramp Drag 

Though the ramp drag is of little interest for itself, it is of 
major importance for calculating the thrust produced under self-propelled 
00110110115 ٠. 

There is no vay of checking the value of the ramp drag as there vas 
in the case of the bollard pull test. Since, however, the ramp drag is a 
contribution to the total model resistance and therefore to Cp, the use 
of %pSU?* as a nondimensionalizing coefficient for the ramp drags might pro- 
vide a basis for their comparison. 

Figures 10 and 11 show the pressure distributions for both models. As 
with the bollard pull test, there are dissimilarities between these curves. 
Again Py is consistently higher on the large model than on the smaller one 
and Py and P3 have exchanged positions relative to Py in the same manner 
5777 the bollard pull test. This flow dissimilarity reinforces the con- 
clusions drawn from the bollard pull test. Namely, that great care must 
be exercised in separating. the forces vhich are acting in the channel so 
that each force is scaled in accordance vith its appropriate scaling lav. 


Figure 12 shovs the ramp drag coefficients, Cup» defined as 


Rp 
b” 1:0 SUZ 


plotted against Froude number. Tables6 and 7 tabulate the original 

data and the values of the intermediate quantities vhich vere calculated 
in obtaining Cen: 
Two facts are apparent. First, there is surprizingly good correla- 


tion considering the basis upon which the values were calculated. Secondly, 


by implication of the good correlation, the value of the force is not 
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Figure 10 


Bare Hull Ramp Drag Pressure Distribution 


as a Function of Position for Model A 
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Figure 1l 


Bare Hull Ramp Drag Pressure Distribution 


as a Function of Position for Model B 
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Figure 12 


Ramp Drag Coefficient, Cap? vs Froude Number, F 
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terribly sensitive to the location of the peaks and valleys. (See 
Figures 10 and 11.) (Later it will become apparent, however, that the 
value is very sensitive to the magnitude of the peaks and valleys.) 

The dissimilarities between these two curves will unfortunately be 
reflected through the rest of the work, since the ramp drag is one con- 
tribution to the thrust calculation. This will be especially true for 
low thrust conditions, i.e., for low air flow rates, because at low air 
flow rates the calculated force which is the difference between the two 
large numbers is more sensitive to errors in those numbers because they 


are close to being equal in magnitude. 
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D. “opelled Hodel Tests 

As explained in the preceding sections on procedures, the self- 
propelled tests are required to provide three additional pieces of in- 
formation vhich are: 

7) props rica speed as a function of the air flow rate. 

(2) The apparent resistance of the hull with air flowing, as meas- 
£— by the force detectors, Tn. 
(3) The pressure distribution on the ramp from which the propulsor's 


thrust can be calculated. 


Results (1) and (2) above are produced simultaneously by simply plotting 
the average measured forces against the air flow rate. 

Figures 13 and 14 present the average value of the measured forces. 
The data is tabulated in Tables 8 and 9 . The zero velocity line does 
not extrapolate to zero by straight line extrapolation. This indicates 
that very low air flow rates produce comparatively higher forces than do 
the higher air flow rates. This is probably the result of a change in 
the mechanism by which the air flow gets converted to thrust as the air 
flow rate increases. It has already been postulated by os D 
that very low air flow rates should be more efficient in the sense that 
the incremental thrust per incremental air flow rate will be high. The 
reasoning is that at low air flow rates the bubbles can mix with the water 
better and therefore the momentum transfer 1s more complete. This work 
seems to confirm his postulate. 

Taking the zero force intercept of the Tg Versus air flow rate curves 


and nondimensionalizing the values to @ and F gives the values plotted in 
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Figure 13 


Self Propelled Test: Apparent Hull Resistance, Tə 


` as a Function of the Air Flow Rate for Model A 
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Figure 14 


Self Propelled Tests Apparent Hull Resistance, Tus 


as a Function of the Air Flow Rate for Model B 
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Table 8 


Self Propelled Tests; 
Apparent Hull Resistance, T 


M 
Model L = 4,5ft 


Air Flow Rate [crn] 15 30 45 60 75 90 


Velocity 0,000 KTS 
T. [ape | 012 085 ,122 .159 195 .237 


Velocity .310 KTS 
m lər) (023 054% „о ,115 is „175 


Velocity .610 KTS 
m, [ies] -.016 „о: „оз .060 6096 „119 


Velocity „901 KTS 
r [be] -.070 =,009 -.031 -.002  .025 .097 


Velocity 1,201 KTS 
T hər) - 119  -.122  -.102  -.090 =,071 -.0/2 


Velocity 1,510 KTS 
Ty | 10£] -.285  -.256 -.204 -.231 -.221 -.192 


Velocity 1,802 KTS 
ту | 102 | - 138 -.827 -.16  -.106  =.393  -.381 
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Table 9 


Self Propelled Tests; 
Apparent Hull Resistance, Ty 





Model L = 6.75ft 
Velocity 0,0 KTS m 
Air Flow Rate [cfh] 148 296 115 593 
HE [1f] 458 —: 1.069 1. 
Velocity 1,006 KTS | ㆍ 
Air Flov Pate[cfh] 153 311 140 589 
ir 1) - 2036 .205 ‚123 ‚681 
Velocity 1,351 KIS 5 | 
Aix Flow Rate [cfr] 148 304 ly? 589 
: E 5 E ar 
7 [uve | ‚305 072 54 3 
Velocity 1,510 KTS — 
Air Flow Rate[cfh] 148 296 149 593 
| RN 0: - 030 ‚166 
m [1er] əl 2 3 
Velocity 1.645 KTS — 
Air Flow Rate [efh] 163 308 160 589 
T. İr) 2.569 ~ =.315 ии ии 
Ы ос ту 1.995 KT»  — 
Air Flow Rate Ffcfh) 148 304 467 50 
-. 8 = 8 -,6 -,519 
n [tef] 98 785 35 





ə 
Figure 15. From this presentation of the data it is apparent that for a 
E... Froude number (velocity) the larger model requires a higher 
value (air flow rate). The plot shows that the self-propelled velocity 
measured progress in a smooth manner and the nondimensional formulation 
allows both lada to be plotted on the same graph. These curves are 
nothing more than plots of the self-propelled velocity versus the self- 
propelled air flow rate. Facts to remember when looking at the plot are: 

(1) That the large model does actually require larger forces at a 
given Froude number because of the tank blockage. 

(2) That the smaller model is in laminar not turbulent flow at the 
low Froude numbers. 

(3) The change in the hull resistance resulting from the air flowing 
has not been separated out of the total forces yet. 

Before continuing to the third set of results derived from the self- 
propelled tests, it should be pointed out that Kostilainen's methods of 
analysis are not very useful for this set of models. Figure 16 shows how 


Cm, varies with Froude number and it is apparent that the correlation is 


TA 
not.very good, even at F = 0 where Kostilainen found good correlation on 
his models. 

One plausible explanation for this is that for the geometry of this 
configuration and the air flow rates used, the nondimensional analysis and 
the assumption used in eliminating some variables.is incorrect. The force 
which is most likely to be causing this is the friction force. (See Sec- 
tions V.B and V.C.) The models used in this work were all tested at lower 


Reynolds' number than Kostilainen's models were and were therefore in a 


flow regime which was more dominated by friction. 
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Figure 15_ 


Froude Number of the Self Propelled Velocity, F, 


as a Function of Nondimensional Air Flow Rate, 
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Figure 16 


Coa vs Froude Number, F 
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The third set of results to be extracted from the self-propelled 
tests is the propulsor thrust, T. Tables 10 through 27 tabulate these 
values and values of other quantities of interest derived from the thrust 
such as Сес , n and E 

With the exception of the fact that an increase in air flow rate is 
normally associated with an increase in the thrust, comparison of the 
thrusts with each other without nondimensionalizing them gives no new 
facts nor does it provide any insight into understanding the underlying 
principles of operation of the propulsor. Section E will take up the 
task of presenting the results in nondimensional form and interpreting 


the meaning of the plots. 
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E. Thrust Coefficient, Air Flow Resistance Coefficient, Efficiency 

The four preceding sections provide all the information needed to 
develop the desired output from this work, i.e., the basis for prediction 
of full scale ship propulsion performance using model data. The theory 
for accomplishing the extrapolation was outlined in Chapter III. 

Kostilainen has suggested that the factor 1/pgB,_h* be used for non- 
dimensionalizing and correlating results of tests on the two-phase pro- 
pulsion system. As shown previously (See Chapter V. D.) this approach 
12 "it give good correlation vhen used to obtain Cpa from T,. Figure 17 
shovs the results of using Poe to nondimensionalize T to yield Cro. 
The correlation is not as good as for Cmp in Figure 16, but the degrada- 
tion may be in part the result of inaccuracies in calculating the thrust, 
T. Elimination of suspicious data points did not significantly improve 
the situation nor did the attempt to perform a form of calibration of the 
thrust calculation based on using correction factors obtained by taking 
the ratios of the measured. forces to the calculated forces in the bollard 
Pal tests. 

It must be concluded that Kostilainen's dimensional analysis is not 
very effective on the geometric form used for this thesis. 

The approach proposed in the theory section is to obtain the hull 
resistance by use of Cp, Cp, and Cy: CH is plotted in Figure 18. Correla- 
tion is very poor and is especially so for air flow rates above the self 
propelled air flow rate. This is because the values of Ce are very sensi- 
tive to value of the total thrusts and some of those values are in error 
as was seen in Figure 9. The sensitivity comes about in the following 


manner. In the calculation procedure for obtaining Cg one of the steps 
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Figuro 17 
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Sto subtract T, from T. At air flow rates above the self-propelled 
point this value should yield a positive number but because T has too 
small a magnitude as a result of the inadequacies of the procedure used 
in its calculation the difference between T and Ta becomes negative. 
Mee, Of course, implies that Cp is negative, since 

E (CTT శతు. 

When Cry is subtracted from the alreddy negative CO E ии. CH the re- 
sultant Cy is very large and negative. 

Below the self propelled air flow rates the magnitude of the Cg 
values will still be too high for any calculation which used T which was 
too small. Fortunately there are three curves in the plot which may be 
reasonably close to the correct value of Cy . These curves are the Cə - 
.146, .209 and .252 curves which correspond to air flow rates whose bol- 
lard pull test integrated force was within 10% of the measured value. 
These three curves cross over each other at a couple of points and one 
value is apparently grossly in error, but they do tend to plot along 
fairly closely. 

The propulsor efficiency for both models is plotted in Figure 19. 
These values suffer from the same maladies as the two preceding quantities. 
Some of the points are obviously in error and can be eliminated from fur- 
ther consideration. Their removal can be confirmed by using the two pre- 
ceding plots to help in the elimination. Unfortunately some points look 
"ad" only on one plot and therefore cannot be totally disregarded. The 
large model's data points plot lower than the L = 4.5 ft. model's data in 


the range of F = .17 and below. This is probably because the calculated 


merees are too low. 
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Once again three curves can probably be believed. They are the 
same ones as in the Cə discussion, i.e., $ - .146, .209 and .252, and 
their credibility is again based on the closeness of the calculated 
forces and measured forces in the corresponding bollard pull tests. 
Based on the selection of those three curves, the efficiency can be 
claimed to be in the range of 20% to 22% and apparently almost independ- 


57 OL air flow rate. 
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VI. DISCUSSION AND CONCLUSIONS 


The conclsuions concerning this work can logically be divided into 
two categories, which are those concerning the experimental techniques and 


equipment and those concerning the results of the investigation. 


A. Equipment and Procedures 

The only equipment which needs improvement is the model itself. Be- 
cause the original model was of the geometry shown in Figure 2 the fol ow- 
on models also had its inadequacies. From the point of view of testing the 
model, the model has too large a beam compared to the breadth of the pro- 
pulsor channel. The extra beam causes added drag, thus reducing the ship 
speed for a given air flow rate and therefore causes a reduction in the 
range of the Froude number in all plots and reduces the value of the maxi- 
mum self-propelled Froude number. 

The air chamber at the bottom of the channel is a waste of ship 
length and should be shortened. A study would be required to determine 
exactly how short that chamber should be and what its real effect on per- 
formance is, but, visual observation indicates that the water flow is be- 
ing restrained from entering the channel at the after end of the chamber 
and this is probably reducing efficiency slightly. 

The fact that the integrated pressures do not equal the bollard pull 
forces indicates that the assumption that the pressure distribution across 
the channel is flat is in error. Visual observation confirms the fact that 
air flow and water flow change as a function of the athwartships position 
17 The channel. In part, this is a result of the fact that the air and 


mater at the higher air flow rates do not mix very well. A new orifice 





ఊం. = 
system might improve the mixing and thereby raise the efficieney Of the 
E opulsor. | 

One possible way of increasing the mixing of the air and water would 
be to increase the air velocity leaving the orifice. Of course this 
would also increase the orifice losses which would tend to reduce the ef- 
ficiency, but a trade-off study could be made to find the optimum air 
flow rate and air velocity condition. 

It should be pointed out that items discussed above are not going to 
result in major changes in the efficiency of the two-phase propulsion 
system. 

The major procedureal inadequacy of this work results from the selec- 
tion of the option of calculating the thrust from the pressure distribu- 
tion vice direct measurement. It was a poor selection for two reasons, 
the first being that the pressure distribution is not flat across the chan- 
nel. In addition to that, the athwartships distribution is probably not 
E e at the bottom of the ship as it is at the stern, therefore, to 
measure the pressure accurately would require a large number of pressure 
transducers and would still probably not be exact or even as accurate as a 
direct force measurement. 

The second reason that measurement of pressures is a poor tool in 
this application is that the pressures are highly oscillatory and produce 
highly oscillatory forces thus requiring averaging and thereby introduce 
errors into the calculations. A significant improvement would be made by 
finding a method for taking longer samples that could yield better average 


values. Placing the propulsor in e flow channel would accomplish this. 
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Therefore it is concluded that: 

(1) A nev geometry be determined for the — 

(2) The ramp be designed in a manner which will allow direct meas- 
urement of the thrust, 

(3) The measurement of the thrust be accomplished in a flov channel. 

These changes will improve the accuracy of the results and shorten 
the length of time required to make the measurements and obtain the data 


needed. 


b. Results 

The results of this investigation are such that few conclusions con- 
cerning the original objectives can be drawn from them. It is obvious that 
many of the points in the plots are in error. The fact that the integrated 
forces do not generally equal the measured forces in the bollard pull tests 
indicates that besides the random errors from mistakes and from the inac- 
curacies of taking averages of oscillating variables, there is also a sys- 
tematic error of considerable magnitude. As noted previously, this sys- 
tematic error probably arises from the erroneous assumption that the pres- 
sure distribution is not a function of the athwartships position. 

Some encouragement can be taken from the fact that when the obviously 
bad results are disregarded, the curves have the same general shape and 
that where they predict nonsense the reasons for the poor quality of re- 
sults can be explained. 

The plot of Cro versus F, Figure 17, gives the best proof that all the 
forces involved in the two-phase propulsion system are not being properly 


scaled by techniques similar to those of Kostilainen. This conclusion is 
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based on the obviously poor correlation between models. 

ine plots of CQ and n versus F, Figures 18 and 19 are too badly con- 
fused with inaccuracies to allow any definite conclusions to be drawn con- 
cerning whether or not they represent the correct way to scale the two- 
phase propulsion system. Both sets of results give some encouragement 
though. Throwing out the obviously bad points, the curves tend to agree 
with physical reasoning. For example, С̧фә displays a behavior which could 
be explained as follows. As the velocity of a vessel increases to higher 
values, the ability of a fixed air flow rate to reduce the hull resistance 
decreases because the increase in water turbulence will increase the prob- 
ability of water penetrating the layer of.air and thereby contacting the 
hul and thus increase the momentum transfer from the hull to the water 
and therefore raise the resistance. This would imply that Co (being nega- 
tive) should move toward positive values with increasing F, and it does. 
Even the "bad" points don't significantly violate this trend. 

It can be concluded that the “¿ and n approach to scaling have not 
been proven to be useless and may be quite good. Better data is needed be- 
fore a stronger declaration can be made. 

Another conclusion vhich can be dravn from this vork is that the ef- 

. ficiency is poor enough so that the two-phase propulsion system will not be 
economically competitive with a displacement ship with a screw. The pcor 
propulsor efficiency will require that extra fuel be provided which in- 
creases operating costs and reduces cargo carrying capabilities. Additional- 
ly, it is apparent from the pictures and the lines drawings that for its 
length it has a relatively small cargo space available, thus increasing the 


initial procurement cost/pound of payload figure. Therefore, on a volume 





US 


limited ship design project, this propulsion system will be of little 
value. Some gains in maintenance costs may be available but the considera- 
tions above will probably disqualify any such gains. 

In summary, it is concluded that: 

(1) The two-phase propulsion system configuration used for this work 
is not very efficient at the air flow rates and velocities measured. 

(2) The scaling law postulated in Chapter III is a potentially useful 
calculational tool and represents an approach to the scaling problem which 
may ultimately be successful. 

(3) Further work should be conducted on the two-phase propulsion 


system in accordance with the suggestions in the next chapter. 
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VII. APPLICATION, ENGINEERING CONSIDERATIONS AND SUGGESTIONS 


There are many other aspects of the two-phase propulsion system 
which must be considered in making a judgement concerning its future use- 
fulness besides the fact that it is not very efficient and will probably 
have a total life cycle cost which is well above that of a conventional 
screw propulsion system. 

Of considerable importance is the problem of backing down. The device 
has no way by itself of creating a backing force. Therefore it would re- 
quire some additional device for that purpose. Kostilainen's self-propel- 
led model solved this problem by providing a ramp on the forward part of 
the vessel to which the air could be ducted and thus provide a backing 
force. Though that solution was satisfactory for a research device it 
would not be practical for commercial cargo-type vessels as it further re- 
duces the space available for cargo. 

Closely coupled to the backing problem is a concern for the vessel"s 
maneuvering capabilities. Most vessels which require any significant capa- 
bility in this respect have rudders which are in the propeller wash. Be- 
se of the highly oscillatory nature of the flow, placement of the rudder 
in the propulsor wash could result in a need for the rudder to withstand 
very high force levels. Placement of rudders on the trailing edge of the 
skegs will undoubtedly be the most practical method of control but will 
probably not result in very good control nara Nn 

The highly oscillatory nature of the force was quite evident at high 
air flow rates. Under those conditions the models experienced considerable 


surging and pitching. If a full-sized ship were forced to undergo such 
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surging and pitching it vould be a very uncomfortable ship to ride and 
nt possibly need structural reinforcement to be able to withstand the 
effects of stress and fatigue. 

Measurements of the pressures on the ramp indicated that for the 
smaller (L = 4.5 ft.) model at high air flow rates the peak pressures 
were commonly as much as twenty times the average. This means that the 
ramp experienced a pounding force which was approximately one-sixth of 
the static head. This situation seems to compare to cavitation which can 
cause severe erosion damage. 

The two-phase propulsion system does have some unique features which 
may ultimately result in its practical application. One is the lack of 
an external screw. A second is that the engine and fan installation would 
be much lighter than that of a conventional power plant. 

One area of application where efficiency is of little importance is 
the use of this device as a lateral thrust unit. It could be used on the 
surface to provide thrusts for keeping ships from putting too much strain 
on mooring lines in high crosswinds by placing a ramp device with camels 
outboard of the ship. 

Along the same lines another possible use of this device might be in 
providing thrust for moving and holding or dislodging objects during shal-- 
low water salvage, repair or construction operations. Ramps could be 
lowered to the bottom and attached to the object to receive the thrust and 
then through the use of compressed air supplied through hoses, the device 
could be operated and controlled. 

It might be possible to use the device on small barges which are towed 


in long lines by a tug. Slight modification of each barge and installation 
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Of one of these devices could reduce the toving pover requirements of the 
uc thus allowing it to tow longer lines of barges and possibly reduce 
the overall cost. 

Continued study of the two-phase propulsion system is recommended 
but the impetus for such work should not arise from an interest in im- 
mediate application to ship propulsion, unless one of the special features 
of the system is anticipated to be particularly useful. There is, however, 
considerable ground for academic interest and for work aimedat establishment 
of fundamental information and also for establishing some of the properties 
of the System which have not been measured. 

One area of study which is quite general and certainly needs some em- 
phasis is the measurement of Forces whose oscillations are of about the 
same magnitude as the average. In the special case of the two-phase pro- 
pulsor, it might be possible to design a force block to attach to a float- 
ing ramp which has one (fore and aft) degree of freedom but is sealed for 
watertight integrity. 

There are quite a number of aspects of this propulsion system which 
have not been looked at yet. Опе such aspect is alternate configurations. 
It would be nice to have a theory which was definitive enough to present 
some suggestions to the experimenter on how to optimize the system through 
geometric variation; however, none exists.' Alternate configurations must 
be approached then from the point of view of applying sound engineering 
principles to the system and establishing what parameters might have the 
greatest effect. Some of these changes are: the ramp angle, the shape of 


the stern, the number and location of orifices and the shape of the ramp. 
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Two aspects of the flow that have not as yet been studied are: what 
happens to the forces when the water flow is in opposition to the bubble 
motion and what is the effect of constraints 하기 lov clearance betveen 
mae bottom of the hull and the bottom of the body of water. 

Though 5 been stated before, it bears repeating that any future 
vork on this propulsion system should not try to measure the thrust by 
pressure integration until direct measurement techniques have been shovn 


to be inadequate. 
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APPEND ZA 


Kostilainen ” has performed a dimensional analysis; a momentum 
balance and an energy balance on a two-phase propulsion system. A sum- 
mary of the analysis is included in this appendix for easy reference. 

It should be noted that the notation used in the following paragraphs is 
the one defined in the main body of this thesis and not Kostilainen's 
original notation. In addition, he deals with a slightly different ge- 
ometry than the one used in this thesis. The main differences are that 
the orifices of models A and B are the full width of the channel, while 
Kostilainen's models had several round orifices evenly distributed across 
the channel, and Kostilainen's models do not have an alr chamber. 

For his dimensional analysis Kostilainen selected as a starting point 


an equation of the form: 


= Flo, Ur > Va» lin. $6» ha Bə H, 0) 


Мр = number of orifices per unit width 
6 = diameter of an orifice 
u = viscosity of the liquid 
g = surface tension of the liquid 
Using the Buckingham 1 theorem he eliminates o, g and h. He then 
presents arguments based on SN results which allow him to 
neglect the dependence of T/B, on o and 1. He then selects 


T 


e = >. 
ogBeh” 


ES 
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Yo 
x " Bohvgh 
and 
İnh 3 ... 
Vgh 


partly based on the fact that these coefficients gave him the best cor- 
relation of data. фу is a dimensionless air flov rate and bar is a 
modified Froude number. 

To obtain a momentum balance he has selected a control volume as seen 
in Figure 20. Using the velocities as shown there (U, is the outflow 
velocity of the gas and U, the velocity of the two-phase flow at the water 
surface) and defining Vi, as the amount of liquid that gets accelerated, 
and B as the angle the ramp makes with the water surface, conservation of 
momentum then gives the ihrust 

B ov. 15 Ur) cos B 
where gas momentum has been neglected. For his energy considerations he 
defines an efficiency, E which is the ratio of the increase of kinetic 


energy of the liquid to the loss of potential energy of the gas. Then 


me 


Лә дк D 
11” xh Ve “027 ÜL”. 


Combining these tvo results he shovs that 


20 Va gh 


= кк CECE 
pU UE 


This could be nondimensionalized. 
Though this has resulted in a rather simple algebraic expression it 
has not made experimentation any easier because ЕТ is not a directly meas- 


urable quantity and U, is at best difficult to measure. 
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APPENDIX B. 


The folloving is a numerical example of the calculation performed on 
each self-propelled velocity-air flow rate combination. (See Tables 18 


and 19.) 


Model: Lh 2.6.75 It. 
Velocity = U = 1.006 KTS 
F= U//eL = .0678 x 1.006 = .1151 


Air Flow Rate = 440 [cfh] (average of several runs) 


Er Tilon - 


L2 VgL 


.001489 x HHO = „655 


Sampson's Rule for calculating the area Under Er eos MID 
Pas P4 and P, yields the following equation. 


F, = force on the ramp parallel to the ship's axis. 


. e 





ir area saa ııı... D 0360 Ж 
ШИ = (.7282)(.2675) = .1948 Ibe 
. 2 

0360 ра BCAPy 2 
S Corr “ —— x ———— 2 .2681 P = (ర í 
B c in H.0 1.5 4 
Bar = Er + Step Corr p ТБТ 
Ilust = T = Total - Ramp Urag 5:5 153 

Look up ramp drag in Table 7. 

Pover in = Va '" pgh =7.0 122 Wale mean ın. 

G G sec 

A : GE 165115 

Power out = T* U= (.3653 1bf)(1.700 ә BT 0 


“The small model was calculated slightly differently. See end of this Appendix. 
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ln is taken from Figure 17. 


“20SU 2C Т - Т = .3653 – .4258 = ㅋ ун 


1 
!49SU 20 


ә 7: 
1... 


2508U “Cpyy = Bare Hull Resistance 
Cpu is taken from Table 3. 
Therefore: 


050206 = – .0577 – .352 = – .4097 lbf 


-.4097 1bf 
ф 2 Tr mas = Oa 
49.07 Ibf 


For the smaller model (L = 4.5722.) chezpre u 5 И 


= у 
vere plotted and a smooth curve was drawn through the points. Simpson"s 
Rule was then cpplied to that curve and was based on a twelve divisions 
vice the four divisions displayed above. The thrusts obtained by making 
the extra plot are within a few percent of the much simpler procedure used 
above and therefore it was not used for the larger model calculations. 

The bollard pull forces and the ramp drag forces were calculated us- 


ing applicable sections of the above procedure. See Tables 4, 5, 6, 


576 7. 





bac 


Flgure 20 


Liquid surfaco at rest 


| 
— 
Momentun Balance for a Tvo-Phaso Propulsor 














Thesis i37 
C75957 Cox 194536 
A systematic study of 


a novel tvo-phase ship 
propulsor. 


thesC75957 


DUDLEY KNOX LIBRARY - 





